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Abstract 

Time evolution of radial wave packets built from the eigenstates of Dirac equation 
for a hydrogenic systems is considered. Radial wave packets are constructed from the 
states of different n quantum number and the same lowest angular momentum. In 
general they exhibit a kind of breathing motion with dispersion and (partial) revivals. 
Calculations show that for some particular preparations of the wave packet one can ob- 
serve interesting effects in spin motion, coming from inherent entanglement of spin and 
orbital degrees of freedom. These effects manifest themselves through some oscillations 
in the mean values of spin operators and through changes of spatial probability density 
carried by upper and lower components of the wave function. It is also shown that the 
characteristic time scale of predicted effects (called Tis) is for radial wave packets much 
smaller than in other cases, reaching values comparable to (or even less than) the time 
scale for the wave packet revival. 
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1 Introduction 



For more than fifteen years large efforts [1-7] liave been made to a detailed understanding of 
the quantum dynamics of wave packets in simple systems like H atoms, hydrogenic atoms as 
well as in simple molecules. These theoretical investigations resulted in a good understanding 
of such subtle interference effects as collapse, revivals, fractional revivals of wave packets 
created in a variety of systems. With development of experimental techniques allowing to 
tailor many different desired initial states many signatures of these phenomena have been 
observed. A rich survey of early studies is given in review papers [8] . 

Radial wave packets (RWP) can be relatively easily excited by short laser pulses [9]. 
Their motion, at the beginning, resembles a classical motion. During later stages of the 
evolution packets undergo dispersion and (partial) revivals as well as fractional revivals 
[10, 11]. The idea of construction of the RWP in the paper is similar to that presented in 
articles cited above. However, we include in our investigation the spin degrees of freedom 
which motion for RWP can manifest itself much earlier than for other kinds of wave packets. 
The most natural framework for our considerations is using relativistic wave packets built 
from the solutions of Dirac equation for hydrogenic systems. Such approach was used by us 
in [12], where we investigated circular wave packets and in [13], where we considered elliptic 
wave packets, both in hydrogenic systems. Results of [12-14] show that contributions given 
by small components of wave functions are negligible small. Therefore for time scales of 
motion which are relevant one can safely use an approximation in which the time evolution is 
calculated with non-relativistic wave functions and relativistic energies. Such approximation 
is simpler for analytical presentation as well as for numerical calculations. We will use this 
approximation throughout the paper. 



2 Construction of wave packet 



Assume that just after creation the radial wave packet has the following form 




(1) 



2 



where \nll) are eigenstates of the non-relativistic hydrogenic system with low angular mo- 
mentum (usually assumed as Z = 1) and m = I. The weight coefficients Wn = (— l)"'c„ = 
(—1)" {2na'^)^^/^ exp [—{n — UavY /^(t"^] are given by a Gauss distribution with mean nav and 
dispersion a. Distributions of that type describe population of Rydberg states excited by 



a short laser pulse. The phase (—1)" is added to obtain the initial localization of the wave 

(a\ 

packet at its external turning point. The spinor determines the initial direction of the 

spin. The WP (0) is an approximation of the full relativistic bispinor whose initial small 
components are set to zero. 



3 Time evolution 



After transformation to the basis \n,l,j,mj) one obtains (with notation j> = I + 1/2 and 
J< = 1- 1/2) 
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In the basis \n,l,j,mj) time evolution of each state is given by an exponential factor 
exp{—iE^it/h) or exp{—iE^^it/h), where E^^ and E~i are energy eigenvalues for j> = I + 1/2 
and j< = / — 1/2, respectively. Precisely 

1 -1/2 



1 + ^ 



{Zaf 



n-j> 



(3) 



Applying time evolution in that basis and transforming back to the \n,l, s,ms) basis one 
obtains the wave packet after time t in the form |\E'r(i^)) = , where the upper \I'i(t) 

Vl^2)/ 

and the lower \E'2(t) component of the spinor are given by 

n 
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^2(t) = 5^ ^/^n I & ( exp {-lE^^tlh) + 21 exp {-lE-^tlK) ) 
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Such wave packet is localized only in the radial coordinate, hence the radial density proba- 
bility of the components 



and P2{r) = r' 



J (in 1^2 
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are convenient quantities illustrating the wave packet motion. The integration over angular 
coordinates leads to formula 
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where Rn{f) denotes the radial part of the wave function {r\nlm). 

Fig. HI illustrates the time evolution of the radial wave packet with rzav = 80, a = 0, 6 = 1, 
what corresponds to initial spin antiparallel to Oz axis. The wave packet exhibits a kind of 
breathing motion, moving towards the center and reassembling itself (approximately) after 
one classical period at the external turning point. 

Periodicity of the motion is well seen with the help of the autocorrelation function A{t) = 
(\l'(t)|\l'(0)). For radial wave packet it reads as 



A{t) =J2wlUa' + b' exp HE+t/h) + b' exp (-^E" 
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The plot of |A(t)P for wave packets with = 80, Z = 92 and two different values of a is 
presented in Fig. |21 

Even for short time evolution presented in fig. Qone can see a transfer of the probability 
density from one component of the spinor to the other already after one classical period. 
It shows that the time scale of the period of spin-orbit motion is for RWP substantially 



smaller than for circular or elliptic WP. The period of spin-orbit motion is determined by 
the splitting of energy levels for which n = ri;^^ is maximally populated in WP and spin 
projections are opposite 
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The result is obtained using lowest order approximation for relativistic energies in a hydro- 
genic system. It is clear that for RWP, whose / = 1, Tis can be even smaller than Trev, 
particularly for large Z. 

A hierarchy of time scales is defined as in [12]. Writing the energy as function of single 
quantum number n, for n = n^-y we define a hierarchy of times 

1 fd''E\ 27ih 
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For = 1 we obtain the classical Kepler time T^, for k = 2 the revival time Trev and so 
on. Fig. El presents the three time scales T^, and Tis for RWP as functions of riav for 
different Z. It is clear that for Z > 60 Tis becomes comparable to Trev and even shorter for 
> 50. Because lifetimes of wave packets with respect to radiative decay are about two 
orders of magnitude larger than Tis (according to [15] Tis/Traddec ~ 0.06) it is much bigger 
chance to observe effects of spin-orbit entanglement for RWP than for any other WP. 

4 Expectation values of spin operators 



Expectation values of spin operators can be easily obtained using equations (jH). For RWP 
they have a simple structure 
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where Un = {E^^ — E^Jh). The terms containing cos (cu^t) and sin (c<j„t) indicate that 
at the beginning one can expect the spin precession followed by the spin collapse implied 



by nonlinear dependence of frequencies Un on n. This behaviour is clearly seen in the 
upper part of fig. ^ where for t/Tis G (0,5) spin vector makes several rotations while for 
t/Tis G (5,20) stays almost constant. At that times the length of spin vector is reduced to 
<T ^ 0.55 which means that the part of spin angular momentum in dynamically transferred 
into orbital motion. Later, for t/Tis G (20,33) spin revives (at half of Tis^ = (2/3)navTis, 
that is t ~ 26.7 Tis for = 80). The high peak of autocorrelation function and the larger 
length of spin vector presented in the lower part of fig. |3] confirm the predicted time of the 
spin revival. 

The spin precession is accompanied by the revivals of spatial probability density. It is 
visible in the autocorrelation function and spin components (fig. |3]) and in details in fig. El 
displaying the spatial probability density for some particular times. 

The inherent entanglement of the spatial and spin degrees of freedom manifested already 
for short times in fig. Q can be illustrated with the help of quantum carpet - space-time 
plots of the WP evolution [16]. Such space-time plot presenting time evolution of pi and p2 
separately is shown in fig. El One sees that if the initial WP has only the \&2 component 
there is a transfer of probability density to the other component and back. This transfer is 
governed precisely by Tig time scale. 

5 Conclusions 

We have discussed the time evolution of RWP in hydrogenic systems using a suitable ap- 
proximation of relativistic approach. The main relativistic effect is the appearance of the 
new time scale due to the spin-orbit coupling. As shown above this time scale can be much 
smaller for the radial WP than for previously discussed cases of circular [12] or elliptic [13] 
WP. This fact implies that in principle experimental observations of some spin-orbit effects 
may become possible with existing techniques. 
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Figure 1: Short time scale evolution of the radial wave packet with n^v = 80, a = 0, 6 = 1 
and a = 2. 
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Figure 3: Log-log plot of time scales (in seconds) Td, T^ev and T\s of radial WP as function 



of for different Z. 
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Figure 4: Time evolution of spin expectation values for RWP with Uav = 80, 1 = 1 and 
a = b (upper part) and the square of the autocorrelation function and the 'length' of the 
spin vector (lower part). 
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Figure 5: Radial probability density p — pi + P2 for RWP with riav = 80, I — 1 and a — b 
for several time instants corresponding to big values of autocorrelation function (revivals). 



Figure 6: Time evolution of radial probability densities pi and p2 for RWP with n^v — 80, 
Z = 1, a = 0, 6 = 1. 
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